Static and dynamical properties of a magnetic impurity in a 
strongly correlated electronic system 



K. A. Hallberg 

Max- Planck- Institut fur Physik komplexer Systeme, Bayreuther Str. 40, Haus 16, 01187 Dresden, 

Germany. 

C. A. Balseiro 

Centro Atomico Bariloche and Institute- Balseiro, 8400 S. C. de Bariloche, Argentina. 

(February 1, 2008) 

Abstract 

Using numerical techniques we study the zero temperature properties of a 
substitutional magnetic impurity in a one-dimensional correlated system de- 
scribed by the Hubbard model with repulsive interaction U. We find that, 
while the static spin correlation maintains a 2kp periodicity, the charge cor- 
relations show Friedel oscillations with a 4fcp periodicity for large values of 
U. The static on-site susceptibility decreases as the interactions in the chain 
are turned on and we interpret this as an increase of the Kondo temperature. 
For large values of U the short range antiferromagnetic correlations dominate 
over the collective Kondo singlet formation. 
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Recently much effort has been devoted to the understanding of the effect of impurities 
in strongly correlated systems However, while the problem of impurities in non- 

interacting systems like the Kondo Anderson and Wolff |J models is well un- 

derstood, very little is known about the effect that the band correlations may have on the 
thermodynamic and ground state properties of impurity systems. Some experiments on 
high-Tc materials fl0|-|l~2f show the effect of impurities interacting with a strongly correlated 



Cu — O plane. For instance, heavy fermion behaviour has been observed in Ndi.gCeo.2Cu02 
below T = 0.3K [10|. This behaviour has been modeled by considering a lattice of 4f Nd 



ions interacting with a system of strongly correlated conduction electrons. The strong anti- 
ferromagnetic correlations in the Cu — O plane is taken into account by breaking the spin 
symmetry [[J . Recent theoretical descriptions of the Kondo effect in a Luttinger liquid have 
been made j|T3| , |14|| finding an algebraic scaling of the Kondo temperature {Tr) with the 



Kondo exchange coupling instead of the exponential relation for the conventional Kondo 
problem in Fermi liquids. 

In this paper we perform a finite-size calculation of a single impurity embedded in a 
correlated electron system. We study the case of a substitutional impurity in a Hubbard 
ring. We consider a magnetic impurity and present results for the ground state properties 
such as charge and spin correlation functions as well as dynamical properties: local density 
of states at the impurity site and frequency dependent spin susceptibility. 

The model Hamiltonian reads: 

H = - *i(4+i,a c i,«T + hc ) + u ^Z n i\ n H + Uin^nii + e ; (n /T + n n ) (1) 

i,cr i^l 

Here Ui is the repulsion on the impurity site I, e/ its diagonal energy, = ti = t' and t. t = t 
for all the other sites. Nevertheless we will consider t' — t — 1 because we find that the 
results are qualitatively the same even if t' -C t. As we are interested in a magnetic impurity 
we take q —t and ei + Ui ^> t. 

We use the Lanczos technique to obtain the ground state of Hamiltonian ([!]) on an A^-site 
chain with N e electrons. In order to get a non-degenerate ground state we consider periodic 



(antiperiodic) boundary conditions for N e /2 odd (even) (closed shell condition). 

In the one-dimensional Hubbard model (in our case it corresponds to U\ = e\ = and 
t' = t) the asymptotic charge and spin correlation functions are given by [|T5[|: 

(n(z)n(O)) = K p /(ttx) 2 + A x cos(2k F x)x- 1 - K " ln~ 3/2 (x) + A 2 cos(Ak F x)x- iKp (2) 
(S(x) • S(0)) = 1/(ttx) 2 + B x cos{2k F x)x- 1 - Kp ln 1/2 (x) (3) 

with interaction-dependent parameters/^, A 2 , B\ and K p . For large U it is expected that 



the 4/cf term of the charge correlation function is dominant [£L6|]. This is confirmed by 



numerical results done in finite systems [17|. As we show below, this behaviour in the 



correlation functions dominates the charge and spin screening of an impurity embedded in 
a one- dimensional system. 

In Fig. [l| we show results for the static spin [(SfSj)) and charge ({niUj)) correlation 
functions for the quarter filled case for different values of U. We see that the interactions in 
the chain enhance the antiferromagnetic correlations and the periodicity remains 2kp for the 
spin correlation. On the other hand, for the charge correlations we find that, as U increases, 
the amplitude of the oscillations decreases and the main periodicity changes from 2kp to 
4k f at U ~ 4t. The behaviour of the system at large U can be understood from the U = oo 
limit were the charge dynamics can be described by non-interacting spinless fermions with 
kp replaced by 2kp. Consequently one expects a contribution proportional to cos(4kp) in 
the density-density correlation functions in this limit. Note that (nf) ~ 1 and {{Sf) 2 ) ~ 1/4 
for all U as corresponds to a magnetic impurity with q <^ —t and q + U\ t. 

We also calculate numerically Jl8| the density of states at the impurity site piioS). In 



Fig. we present results for the doped case where the delta peaks have been broadened by 
a Lorentzian only for visualizing purposes. Important peaks exist close to ei and ei + U\ (the 
latter ones are out of the range of the figure). The lower and upper Hubbard bands can 
clearly be recognized for U = 6t (Fig. 2c). A large density of states near the Fermi energy 
can be seen for all U values. This enhancement is due to the scattering potential introduced 
by the impurity. Due to the finiteness of the system, we can see two enhanced peaks, one 



beneath and one above the Fermi surface but these should become one resonant peak in the 
infinite system. In the Anderson model (in the spin fluctuating limit) this corresponds to 
the Kondo resonant peak and its width is proportional to ksTx- Because we have a finite 
system it is impossible to extract from these figures the width of the resonant peak. However, 
we observe that as U increases, the most prominent peaks separate somewhat; this could 
suggest that the resonant peak gets wider in the thermodynamic limit, implying an increase 
of the Kondo temperature. For the half filled case the resonant peak at U = disappears 
at finite U due to the formation of the Mott-Hubbard gap. The results presented above are 
consistent with the dependence of the static spin susceptibility with the interactions of the 
band presented below. 

The local frequency and temperature-dependent spin susceptibility is given by the fol- 
lowing expression: 

AlK ' W-oo \2k B T) uo-uo' + i5 V 1 

where (( )) R is the retarded Green function and x'i is i n units of (gyUg/2) 2 . 
The static and zero temperature limit is: 
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xj<* = 0) = -2Re MSt u + Eo _ H + iS SrM ^ = 2E™Ir (5) 
where |^o) is the ground state of Hamiltonian ([!]) with energy E and \u) and E u are the 
excited states and their energies, respectively. This static susceptibility describes essentially 
the inverse of the difference in energy between the first excited and the ground state {AE). 
This can be seen by looking at the imaginary part of the response function: 

X '{(u) = -ilm(^\S+ 1 — — STIVb) (6) 

7T UO + Eq — H + ZO 

This is shown in Fig. 3 where we can see that the first excited state has much more weight 
than the other states of higher energy and is clearly distinguished from them. In the half 
filled case, instead, the spin excitation spectrum is more complicated but still a clear first 
excitation peak can be distinguished for small U. We have also checked the spin of the 
ground state and first excited state, finding they are a singlet and triplet respectively. 



For U = 0, i.e. no correlations in the band, Hamiltonian ([!]) corresponds to the Wolff 
model in ID |J. In this case it is known that the properties of the system are very similar to 
those of the Anderson model in the spin fluctuating limit ||19|| . In particular, as calculated 



in J2U] and checked numerically by us, x'i increases with the impurity interaction Ui. 

In Fig. |3] we plot the susceptibility as a function of U for several fillings. The case N e = 2 
(not shown) is special because one electron is localized at the impurity site so there are no 
many-body effects on the band. Provided q <C —t and e\ + Ui t (as is our case), a simple 
second-order perturbation calculation shows that AE always decreases with U. When more 
electrons are added to the system and the correlations start playing an important role, we 
find that x'i has a non-monotonic behaviour as a function of U. For small U the susceptibility 
decreases, it reaches a minimum and then starts to increase. This behaviour can also be 
seen from where the first peak goes towards higher energies for small U but then tends 
to zero (very high x'i) f° r large U (see Fig. |). The susceptibility for the half filled case, 
instead, always increases with U. The inset of Fig. |] shows the susceptibility as a function 
of U for a quarter filled band for N = 8 and N — 12. Here we see that the lowering of x'i f° r 
low U is more pronounced for larger systems suggesting that it will be an important feature 
in the thermodynamic limit. 

In conclusion, we have calculated numerically some properties of dilute magnetic im- 
purities embedded in a one-dimensional correlated electronic system. Although we have 
performed our calculations in a finite system where the low energy characteristics may be 
missed, this method gives an insight on the properties of the system at zero temperature. 

The static correlation functions reflect the non-Fermi liquid behaviour of correlated elec- 
trons in one dimension: While the periodicity of the spin-spin correlation remains 2k f, 
the Friedel oscillations in the density-density correlation have a 4k f periodicity at large U. 
This could probably be measured experimentally with X-ray scattering in doped quasi-one- 
dimensional systems with dilute impurities. 

We also find a large local density of states near the Fermi energy due to the impurity 
scattering potential. The separation of the most important peaks may indicate that the 



width of the resonant peak in the thermodynamic limit increases. 

The susceptibility, which measures essentially the inverse of the energy difference between 
singlet ground state and first excited triplet state (interpreted here as proportional to Tk), 
shows a non-monotonic behaviour for the doped case. For U = a bound singlet state 



is formed around the impurity In this non-interacting limit only a small fraction of 

electrons (those within ksTx of the Fermi surface) are involved in the extended singlet state. 
When the interactions in the band are turned on, the effective mass of the carriers increases, 
enlarging the density of states near the Fermi energy. As more electrons are involved in 
the screening of the impurity, the spin compensating cloud becomes more localized and as a 
consequence, T K is larger i.e. x[ decreases. For large U this picture holds no longer because 
the strong antiferromagnetic correlations inhibit the formation of a collective singlet state 
around the impurity and the susceptibility increases monotonously. In the half filled case 
a gap opens for finite U and there is no formation of a collective singlet state around the 
impurity site and the susceptibility also increases. In this case one expects the susceptibility 
to diverge at large U (when Ui — > oo) because conformal field theory studies ]2T] show that 
the fixed point of a substitutional spin with antiferromagnetic coupling in a Heisenberg chain 
is a 'healed' Heisenberg chain (with J = 2t 2 /U). 

For the highly doped case, our system is in the Luttinger liquid regime and our results 
for Tk at small U are qualitatively similar to those of Ref. (small and positive gi in 
their case and small Kondo interaction parameter J). We find that in the Hubbard model 
the antiferromagnetic correlations play an important role and compete with the collective 
Kondo singlet formation, dominating at large values of U. 

It is interesting to study whether in higher dimensions the features encountered here, 
i.e. the the non-monotonous behaviour of the Kondo temperature as a consequence of a 
competing interaction between Kondo singlet formation and antiferromagnetic interactions, 
persist. Of course, other methods than exact diagonalization techniques are required. The 
doubling of the periodicity of the Friedel oscilations in the density-density correlation func- 
tion is, however, a characteristic of Luttinger liquids. 
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FIGURES 

FIG. 1. Static spin-spin (a) and density-density (b) correlation functions vs. the distance to 
the impurity site for N = 12, Ne = 6 (k F = tt/4), U x = lOOt ,ej = -10* and U = (circles), J7 = t 
(squares), U = 2t (diamonds), U = 4t (triangles) and U = Wt (plus). 

FIG. 2. Density of states at the impurity site I for N = 10, iVe = 6, U\ = WOt and e/ = —Wt 
for U = (a), 2t (b) and Gt (c). Full (dotted) lines: occupied (empty) states. 

FIG. 3. xf( w ) f° r the same set of parameters as in Fig. 1 but for U = (full line), U = 4t 
(dotted line) and U = 30t (dashed line) 

FIG. 4. Local static susceptibility x'i as a function of U for Ui = lOOt, e/ = — Wt, N = 10 and 
different fillings: iVe = 4 (full line), 6 (dashed line) and 10 (dashed-dotted). The inset shows x'i as 
a function of U for a quarter filled band for N = 8 (lower curve) and N = 12 (upper curve), for 
[// = 20£ and e t = - Wt. 
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FIG. 1a 




FIG. 1b 
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FIG. 2b 
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FIG. 2c 
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